Atomistic quantum transport modeling of metal-graphene nanoribbon heterojunctions 
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We calculate quantum transport for metal-graphene nanoribbon heterojunctions within the atom- 
istic self-consistent Schrodinger/Poisson scheme. Attention is paid on both the chemical aspects of 
the interface bonding as well the one-dimensional electrostatics along the ribbon length. Band- 
bending and doping effects strongly influence the transport properties, giving rise to conductance 
asymmetries and a selective suppression of the subband formation. Junction electrostatics and p- 
type characteristics drive the conduction mechanism in the case of high work function Au, Pd and 
Pt electrodes, while contact resistance becomes dominant in the case of Al. 
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Integration of graphene-based nanostructures in elec- 
tronics, sensors and environmental applications makes 
necessary a clear understanding of the interaction be- 
tween graphene and metallic surfacea^"^. Interface bond- 
ing and electrostatics can play a crucial role in the 
transport characteristics of these systems since the low- 
dimensionality and high carrier mobility of the channel 
materiaP can enhance the role of the metallic contact 
with respect to the traditional complementary metal- 
oxide semiconductor technology. In this sense it can be 
argued that the main source of resistivity in graphene- 
based devices should derive from the interaction with 
the metallic electrodes. Characteristics of such inter- 
action for two-dimensio nal g raphene have b een identi- 
fied both experiment ally2EEl anc j theoreticalljEEl, where 
charge transfer, doping-related phenomena and near- 
interface potential fluctuations have been reported. How- 
ever, as patterning and lithographic techniques advance 
towards one-dimensional (ID) confinement in order to 
engineer the necessary bandgaps for digital applications, 
a particular ID electrostatic response can be expected 
that should strongly differentiate device characteristics 
with respect to the two-dimensional case^. Under this 
perspective we study metal-graphene nanoribbon (GNR) 
heterostructures within self-consistent quantum trans- 
port simulations on the basis of: a) an atomistic de- 
scription on both the active device part and the metal- 
lic electrode that respects the interface chemical bond- 
ing, b) a proper treatment of the junction electrostat- 
ics and c) depletion region length-scales. Results show 
that band-bending and doping effects can significantly 
alter the ideal transport characteristics of GNRs giving 
rise to asymmetries in the conductance and a selective 
suppression of the ID subband formation. Moreover, 
electrode-dependent scattering processes can block con- 
duction channels in particular cases. Similar to carbon 
nanotubes (CNTs), we find that there are long-range 
depletion tails in the charge distribution^ that vary on 
the basis of the conductive character of the respective 

gnrP. 
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FIG. 1. Configuration scheme of the simulated graphene 
nanoribbon systems, where a GNR is end-contacted by a 
three-dimensional semi-infinite metallic electrode at the left 
side, whereas ideally contacted at the right side. 



We consider hydrogen-terminated armchair and zigzag 
graphene nanoribbons (aGNRs and zGNRs respectively) 
and use the terminology of Ref. II to categorize them 
on the basis of the dimer lines N a (zigzag chains N z ) 
along the ribbon width. Fig. [I] shows the two-terminal 
geometry used throughout this study, where GNRs with 
channel lengths L c h R3 YJnm are end-contacted at the 
left side by the (111) surfaces of three-dimensional semi- 
infinite electrodes (Au, Pd, Pt and Al). The right elec- 
trode is an ideal ohmic contaclP^, i.e. a GNR with the 
same dimer lines (zigzag chains) as the device part. In 
the case of semiconducting aGNRs this geometry cor- 
responds to a Schottky junction. We employ a self- 
consistent Schrodinger/Poisson scheme for the calcula- 
tion of transport and electrostatics. Quantum transport 
is computed within the non equilibrium Green's function 
formalism (NEGF) coupled to the standard Landauer- 
Buttiker approactH the single particle retarded Green's 
function matrix reads Q = [ES — H — — E^] -1 , 
where E is the energy, H (S) is the device Hamiltonian 
(overlap) matrix and ^l,r are self-energies that account 
for the effect of scattering by the contacts (S = Tg s r', 
where g s is the surface Green function specific to the 
contact type and r is the Hamiltonian relative to the 
interaction between the device and the contact). From 
the total transmission probability T — Trace [T^G^rQ'] 
(where = i^L,R — R ]) conductance can be cal- 
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culated as G — (2e 2 /h)T. The device spectral func- 
tion is the anti-hermitian part of the Green matrix 
A = i(G — G'), from which the local density of states 
(LDOS) at energy E and position r a can be defined 
as: LDOS(r a ,E) = f R3 Trace[AS/(2n)]6{r ~ r a )dr, 
where 5 is the Delta function and r a shows the positions 
of the atomic sites. Hamiltonian and overlap matrices 
are written within a first-principles-base d par ameterized 
model using the extended Hiickel theory-2HSl and a non- 
orthogonal double-C Slater-type basis that fits the band- 
structure of bulk graphene^ and fee metals^ from den- 
sity functional theory calculations. Metal surface Green 
functions for the evaluation of the respective self-energies 
are calculated for the three-dimensional semi-infinite con- 
tact with a back-and-forth real to fc-space Fourier trans- 
form exploiting lattice periodicity^. Charging effects 
are introduced in the formalism with the inclusion of 
a self-consistent potential U sc {pf) that is a functional 
of the device density matrix and is added to the bare 
device Hamiltonian. Within the self-consistent proce- 
dure, mobile charges pf deriving from the NEGF are 
passed to a three-dimensional numerical Poisson solver 
V 2 C/ SC = —pf/e, considering the device part embed- 
ded in SiO'p^. A Dirichlet boundary condition is set 
in the metal-GNR interface of the Poisson box with a 
value U'fJ 1 = 4> rn — <j) gr , where (f> m , (f> gr are the ex- 
perimentally measured work functions for (111) metal- 
lic surfaces and graphene^. Null Neumann boundary 
conditions are set for the other five faces of the Poisson 
simulation box. Self-consistency is enhanced by a pre- 
dictor/corrector Newton- Rapson algorithm^ while opti- 
mized matrix manipulation techniques^ have been im- 
plemented throughout the numerical code. Fermi-Dirac 
statistics are introduced for room temperatures (300K). 

Fig. [2] shows a real-space representation of the band 
formation along the ribbon lengths within total/local 
density of states spectra for a semiconducting N a =16 
aGNR and a semimetallic AT =14 aGNR. In the case 
of the 16 aGNR contacted with the high work func- 
tion Au electrode (Fig. [2ja)) the equilibrium Fermi level 
alignment for the two parts of the heterostructure gives 
rise to significant upwards band-bending phenomena near 
the metal-aGNR interface due to the higher work func- 
tion of the metal with respect to the GNR. However, 
band-bending is not rigid for both conduction and va- 
lence bands as a result of a complex interference mecha- 
nism: the LDOS distribution clearly shows the presence 
of wavelike quantum interference patterns due to the re- 
flection of the incident electron wave by the non-ideal 
contact^. Near the interface such patterns tend to turn 
upwards for the conduction band and downwards for the 
valence band and respond differently in the presence of 
the electric field induced by the barrier. Hence, conduc- 
tion band shifts smoothly while valence band shows local- 
ization patterns in the LDOS distribution. Such patterns 
become discrete localized states with a few-nm spatial 
breadth in the energy region where the bended valence 
band is triangularly-like confined inside the bandgap. In 
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FIG. 2. Density of States (left) and real-space bands along 
the GNR length (right) for (a) a 7V a = 16 aGNR contacted 
with Au, (b) a 7V a = 16 aGNR contacted with Al and (c) a 
iV a =14 aGNR contacted with Pt. AE F in (c) denotes the 
difference between the Fermi levels of the metal-contacted and 
the respective ideal aGNR. 



addition, metal-induced gap states (MIGS), i.e. tails 
of the metallic wavefunctions decaying very fast in the 
semiconducting gap, form throughout the interface (vis- 
ible as a brighter left-border line for all energies in the 
LDOS representation of Fig. [2ja)). It can be therefore 
argued that the interface between a GNR and a metal- 
lic contact is ruled by complex band-bending, interfer- 
ence and localization phenomena whose influence in the 
conduction mechanism will be discussed in the follow- 
ing. When the same aGNR is contacted by a low-work 
function Al electrode (Fig. [2Tb) ) the bands bend down- 
wards {4>ai — <j>GNR < here), whereas qualitatively sim- 
ilar behaviors as before (interference patterns, localized 
gap states, MIGS) can be observed. In both cases the 
Fermi level remains within the bandgap although loos- 
ing the midgap position of the respective ideal aGNR. In 
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the case of a semimetallic 14 aGNR contacted with Pt 
(Fig. [2jc)) the main issue arising from the interaction 
between the two structures is a p-type doping effect due 
to the presence of the high work function metal (see AEp 
in Fig. [5Jc) for the difference between the Fermi levels 
of the metal-contacted and the respective ideal aGNR). 
Hole carrier injection has been obtained for all high work 
function metals on metallic GNRs in this study while a 
less pronounced electron doping effect has been observed 
in the case of Al. Band-bending is also evident here from 
the first 7r — 7r* bands and onwards, however the pres- 
ence of the electrostatic potential does not seem to affect 
the states that lie inside the first ir — ir* plateau (e.g. see 
the GNR-long flat line that corresponds to the secondary 
meV bandgap of the 14 aGNR at the ideal structure's 
Fermi level in Fig. |2jc)). 

Characteristic ID junction electrostatics are present 
in the metal-aGNR case. Fig. [3] shows potential profiles 
along ribbon lengths for the previously shown iV a =14 
and 16 aGNRs contacted by all available metals in the 
study. The main aspect of the electrostatic potential for 
the semimetallic aGNR is a steep potential drop near 
the contact interface that decays after few nm to a non- 
zero flat value. This finite potential value denotes the 
presence of carrier accumulation throughout the GNR 
length (holes for Au, Pd and Pt and electrons for ^4i). 
In the case of the semiconducting 16 aGNR the Schot- 
tky junction behaves qualitatively different. Screening 
is smoother and charges tend to vanish away from the 
metal contact, however also in this case long-range deple- 
tion tails in the charge distribution have been obtained, 
in accordance with previous studies on CNT junctions^. 
In this sense an accurate estimation of depletion length 
scales becomes difficult in these systems and "breaks" 
the traditional metal-semiconductor scheme, giving rise 
to novel ID device design possibilities. The categoriza- 
tion of metal-GNR electrostatics on the basis of the con- 
ductive character of the respective GNR has been also 
encountered in the metal-CNT case^l. It can be argued 
that as the width of semiconducting GNRs grows and the 
respective bandgaps decrease^ we can expect an electro- 
static response that smoothly shifts from Fig. |3jb) to 
Fig. §a). 



FIG. 4. Conductance as a function of energy for a N a ~U 
aGNR (left column) and a JV a =16 aGNR (right column) con- 
tacted with: (a-b) Au, (c-d) Pd, (e-f) Pt, (g-h) Al. Dashed 
lines show ideal conductances for the respective aGNRs. 



Fig. [4] shows the influence of chemical bonding and 
electrostatics in the conduction mechanism of the stud- 
ied systems. High work function Au, Pd and Pt metals 
give rise to qualitatively similar transport characteristics 
that originate from the electrostatic aspect of the het- 
erojunctions. Namely, p-type conduction characteristics 
have been obtained for the 14 aGNR and low Schottky 
barriers with respect to the valence band (of the order of 
0.2-0.3 eV) for the 16 aGNR. Fermi level to conduction 
band distances increase for the semiconducting ribbon 
with respect to the ideal case, arriving at Ec — Ep ~ leV 
for Pd and Pt. In all cases conduction band charge flow 
is strongly suppressed, giving rise to a selective loss of 
the quantization steps that are typical of the ID sub- 
bands in GNR structures. This behavior is related with 
the smooth bending of the conduction band that creates 
a state- free zone near the interface (see Fig. [2]ja)). The 
combination of p-type characteristics and conductance 
suppression due to band-bending gives an asymmetric 
form to the overall conductance distribution (as similarly 
calculated also for CNTs ?1 ). In terms of chemical bond- 
ing only Au seems transparent near the Fermi level with 
the conductance arriving at the IGq plateau of the ideal 
case, whereas Pd and Pt demonstrate a slightly lower 
transparency. On the other hand, valence band trans- 
parency above the first conductance plateau is enhanced 
for Pd and Pt, which show a smaller extent of conduc- 
tance fluctuations with respect to Au, making them more 
appropriate for high bias electrical measurements. A 
careful comparison between group 10 transition metals 
Pd and Pt shows that nonetheless the similarities deriv- 
ing from their electronic structures, Pd shows a slightly 
better conductance response in the quantization steps of 
the valence band. The case of low work function Al 
electrode is distinct, since despite the contact-induced 
n-type doping (for the 14 aGNR) and quasi-ambipolar 
Schottky behavior (for the 16 aGNR), the dominant as- 
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FIG. 5. Electrostatic potential profile U 3C as a function of 
the channel length L c h (a) and conductance as a function 
of energy (b) for a iVz = 10 zGNR contacted with Au. The 
dashed line shows the ideal conductance of the zGNR. 

pect that characterizes conduction is the strong scat- 
tering by the contacts. Here contact resistance consti- 
tutes the main factor of conductance suppression with 
respect to the ideal case, with quasi-blocked conduction 
channels and overall conductance degradation through- 
out the energy spectrum. It is therefore clear that the 
electrostatics and chemical bonding act complementary 
in metal-graphene nanostructures and a categorization of 
the metallic contacts on the basis of their transparency to 
graphene should incorporate a best compromise between 
these two aspects. Finally it should be noted that local- 
ized gap states that form near the metal-GNR interface 
(see Fig. [2]Ja),(b)) do not contribute to the transport 
process. 

Junctions between metals and zGNRs preserve sim- 
ilar qualitative characteristics with respect to aGNRs. 
However, the presence of the edge states in the channel 
material and the accompanying large DOS near the 
Fermi level of these systems strongly enhances the role 
of localized electron-electron interactions. Hence, con- 
trary to aGNRs, the electrostatic response is not uniform 
throughout the zGNR width and gives rise to a faster po- 
tential screening near the borders than in the center of 
the zGNR (Fig. [5|a)). Moreover, the reduced area of in- 



terface overlap between metallic and edge wavefunctions 
further hinders the transparent transmission of electrons 
in these systems (see the lower conductance with respect 
to the ideal case in Fig. [5](b)) . It should be noted though 
that by the suppression of the edge state (e.g. due to cor- 
rugation from nanolithographic processes) , transport and 
electrostatic properties are expected to converge towards 
the aGNR case. 

To conclude, this study has addressed the problem 
of metal-GNR heterojunctions within an atomistic ap- 
proach that deals with both the electrostatics as well 
as the chemical aspects of the interface. Results have 
shown that band-bending, doping and bonding charac- 
teristics of this interaction can non-trivially influence 
the conduction mechanism, giving rise to conductance 
asymmetries, Schottky barriers and suppression of ideal 
transport properties. This study also implies that the 
electrostatics and the chemical bonding aspects can act 
complementary for the determination of contact trans- 
parency in graphene. GNRs, as ID sp 2 carbon al- 
lotropc systems share a lot of common properties with 
CNTs. Within a certain qualitative framework, this 
work argues that theoretical/experimental knowledge ob- 
tained for metal-CNT heterojunctions can be also valid 
in the case of GNRs. It is therefore crucial to under- 
stand the pros and cons of the two systems in terms 
of fabrication/growth/patterning methods and electri- 
cal/mcchanical/optical characteristics in order to distin- 
guish the ideal candidate for post-Si nanoelectronic ap- 
plications. 
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